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Abstract. The square tight-binding model in a magnetic field leads to the almost- 
Mathieu operator which, for rational fields, reduces to a q x q matrix depending on the 
components //, v of the wave vector in the magnetic Brillouin zone. We calculate 
the corresponding Green's function without explicit knowledge of eigenvalues and 
eigenfunctions and obtain analytical expressions for the diagonal and the first off- 
diagonal elements; the results which are consistent with the zero magnetic field case 
can be used to calculate several quantities of physical interest (e. g. the density of 
states over the entire spectrum, impurity levels in a magnetic field). 
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1. Introduction 



Different approaches to the problem of an electron in a two dimensional periodic potential 
with an applied perpendicular magnetic field (see e.g. M or 0) allow the reduction of 
the two dimensional eigenvalue equation to a one dimensional difference equation of the 
form 

9m+i + 9m-i + 27 cos(27ram - v)g m = eg m , (1) 

known as almost-Mathieu or Azbel-Harper equation, where e is the energy of the 
electron scaled in units of the bandwidth at zero magnetic field. 
Equation ([J) corresponds to an effective Hamiltonian 

ii e ff = 2cos(]3 ;c ) + 27cos(27raa; — u). 

It has been studied extensively ever since Hofstadter's pioneering work [|TJ, both for 
the fascinating cantor-set properties and selfsimilarities of its eigenvalue spectrum 
spec = {e|g m non diverging} (see e.g. @, @, H and references therein) and for its 
significance for periodically structured 2d-electron systems (e.g. || and 0). 

In the following we restrict ourselves to the symmetric case 7 = 1 corresponding to 
square symmetry of the 2d periodic potential. 

The phase v e [0, 27r] in (|1|) represents the component of the wave vector in y- 
direction, a is the number of magnetic flux quanta per unit cell or — depending on the 
chosen ansatz — the number of unit cells per flux quantum. 

As in most other work on this problem, we restrict a to the dense set of rational 
values 

V 



a 



q 



With this restriction ([J) becomes periodic, i.e it remains invariant under the substitution 
m 1— > m + q. We therefore may assume m G {0, 1, . . . , q — 1}. According to Floquet's 
theorem ([I]) then has at least one solution of the form 

g m+q = e^g m . (2) 

The phase fi can be identified with the wave vector in x-direction. With (0) the 
solution of the almost-Mathieu equation ([IJ) can be reduced to the eigenvalue problem 

HV> = eip (3) 

of the g-dimensional matrix corresponding to the almost-Mathieu operator H for 
rational a: 
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where Cj := 2cos(27rja — v). 



As has first been shown in ||, the eigenvalues of (Q) can be obtained as the roots 
e of the equation 




Here -P(A) is defined as the polynomial of degree q in A for any A G C: 




Evidently for each fixed pair of values (/z, v) (i.e. one point in the magnetic Brillouin 
zone) (fl) yields q real eigenvalues. Varying v G [— tt, tt] and ji G [— 7r, 7r], the eigenvalues 
are broadened into q magnetic subbands e(/z, v) each eigenvalue e is a function of 
the momenta v and fi and of the subband-index k, with 1 < k < q. Except for some 
special cases where two subbands share one common point, e.g. for even q at e = 0, 
these subbands are disjoint. The graph of this subband structure in dependence on a 
is known as Hofstadter's Butterfly 

According to ([|), the eigenvectors ip^ in ([3D are vectors in C q with the elements 
gi,...,gg. Each is a solution of flU) for e = e/-(//, i/) it depends on z/, /i, on the 
subband-index and on its lattice site j. 

In the next section we are going to develop a method to evaluate the resolvent or 
Green's function G of H without explicit knowledge of the eigenvalues and eigenvectors 
of H. This approach allows us to calculate an analytical expression for the diagonal 
elements of the resolvent in section |3|, results that go beyond those obtained by Ueta in 
0] In section [| we show how to calculate the off-diagonal elements of G. A conclusion 
is given in section 

2. The Green's function 

The Green's function G(A,/i, u) corresponding to H(yU, u) is defined as 



for all A that are not in the spectrum of H 

If the eigenvalues e^ifi, v) and the eigenvectors ipkit^, v) of H were already known, 
G(/io, vq) could be evaluated through 



For most purposes (e.g. the calulation of impurity states) it is not sufficient to 
know G at only one point (fiQ, u ) in the magnetic Brillouin zone; instead the Green's 
function of the original 2d-tight-binding Hamiltonian, i. e. over the whole spectrum of 



G(A, //,*/) = (A -H(//,*/)) 



-i 




(7) 
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H, with fi G [— 7r, 7r] and z/ G [— 7r, 7r] is needed. This is equivalent to a summation over 
the entire spectrum {efc(/i, i^)}, i.e. by integrating over /i and v: 



47r2 t^i i \-e k (n,v) 

— 7T — 7T 

Alternatively the (g x g) matrix G(/i, z/) can be evaluated applying Cramer's rule 
for the inverse of a matrix: 

(G(\,» ,Vo)) m>n = ((A-H)- 1 )^ 

A 



det(A-H) 

= Am < n . (9) 

-P(A) — 2 cos(gz/ ) — 2 cos(gyUo) 

For fixed values /i , ^o- ^m,n in © represents the classical adjoint of A — H i. e. the 
determinant of the matrix that we get when we delete the m-th row and the n-th 
column of A — H, multiplied by (— l) m+ra . 

As in equation (|) we get the Green's function of the original Hamiltonian, that is 
for all possible values of fi and v through integration: 

7T 7T 

Iff A 



G ( A )W = 4^ J J P(A)-2cosM-2cosM d " d/i - 



In this representation G depends only on H itself and not explicitly on the eigenvalues 
or eigenvectors of H. 

The Green's function (|T0|) has some really surprising features: 

(i) An inspection of @ shows that 

G(A)Gl M , (11) 

for all A G R, A G" spec(H), because the imaginary part of A TO>n , if nonzero, has 



always a factor sin(gyu), so it vanishes under the integral (|10 

(ii) Due to the hermiticity of H and (|TT|): 

(G(A)) m> „ = (G(A)) n>m (12) 

for A G R. 

(iii) The elements in every (off-) diagonal are equal to each other: 

(G(A)) n , n+j = (G(A)) lil+j (13) 

for all 1 < n < q and < j < q — n. The proof of this property uses that 

A n>n +j(v) = A lt i +j (i> - 2nan), 

which is easy to accept but rather lengthy to prove ||10|| ; with this result the 
substitution v \— > v' :— v — 2nan leaves the integral (|l0| ) invariant. 
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3. The diagonal elements 

Using equations (|T0| ) to (|13|) we are able to calculate the diagonal elements of G. At 
first we take advantage of 

H 3=1 

From 



ddet(A-H(/i,z/)) d P(A) - 2 cos(g/i) - 2 cos(» dP(A) 



dA dA dA ^ * 3J 



the elements of (|10|) with m = n can now be transformed to 



:G(A)) - = 4^^T / / P(A)-2cosM-2cosM d ^- (14) 



The double integral in ( |T4] ) is well-known, e.g. ||TTJ, and (G(A)) evaluates to 



2 1 dP(A) / 4 \ 

^ G(A)) - = VqPjX) — K [pJXj) 



(15) 



for all A not in the spectrum of H, hence, according to (||), for |P(A)| > 4. K denotes 
the complete elliptic integral of the first kind, cf. ||12|| . 

In order to calculate the density of states (and other useful quantities) we need the 
well-known extension of G to complex A and we define 

G ± (A) = limG(A±is), 

s— »0 

which is relevant for A G R in the continuous parts of spec(H). 

The diagonal elements of G ± (A) are calculated with the analytic continuation of 
the elliptic integral [11]. For the real and the imaginary part we get 

R(G±(A)) =sigB(P(A))^^KOT 

V v jj n!n & v v )) 2nq dA V 4 / 

V v jj n>n ^ 2vrg dA V 4 / 

The density of states p is given by the imaginary part of G + : 

This is identical to the result of Wannier, Obermair and Ray 0. 
The diagonal elements of G and G + are shown in figure [I|. 
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Figure 1. The diagonal elements of G(A) and G + (A) versus A: real (- 
imaginary part ( ) for (a) a = | and (b) a = |. 



4. The Off-Diagonal Elements 

The off-diagonal elements of the Green's Function are required for many purposes, e.g. 
for the calculation of impurity states [[HJ. In the case of a vanishing magnetic field, 
where -P(A) = A, it is possible to find a recurrence relation for all elements of the 
Green's Function [[H|]. In presence of a magnetic field we were not able to find a general 
recurrence. Only the elements of the first off-diagonal (G(A)) nri±1 may be evaluated in 
a similar way: 
Using that 

(A-H)G(A,/i, v) = 1 

and 

G(A,/i, i/)(A-H) = 1 
for every /i and z/, we get after integration 

(G(A)) mj + (G(A)) i _ 1J + (G(A)). . +1 + (G(A)) M ._ 1 = 



7T 7T 



' 1 r (2A - Ci - Cj) — rr ^ -— d/xdi/-25i 



4vr 2 J J v "" ^ P(A) -2cos(g/i) -2cos(g^)' 

— 7T — 7T 

With (|12"D we get an expression for the elements of the first off-diagonal: 



(A - Ci) — -r ^ -— d//dv - 1 I . (17) 

v tJ P(A) -2cos(g//) -2cos(gz/) P 1 




The integral in (|TT| ) yields [|T^ 



(G(A)) 1)4+1 = \ (A (G(A)) m - 1 
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All the other elements of the Green's Function have to be calculated numerically 
by integrating fllQ|) , the result for a — | is shown in figure 0. 




Figure 2. All different elements of G for a = -|: (G(A)) n n+l versus A with I = (a), 
I = 1 (b), I = 2 (c), I = 3 (d) and I = 4 (e). 



For some purposes it ist necessary to extend the Green's Function over more than 
one magnetic period in x-direction (\m — n\ > q). For this case a continuation of G 
can by constructed by applying equation (|j) to equation (0). The elements of this 
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continuation are given by 




-P(A) — 2 cos(qu) — 2 cos(gz/) 



dz/d/i 



where /i is the integer part of (m — n)/q and n' the remainder of n/q. 

For each fixed n' G [0, g — 1] we see in figure ^ with a = | as an example, that 
G(A, m, + /eg) decreases exponentially with growing \m — n| (see figure |^). 
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Figure 3. The extended Green's Function (G(A)) TOJJ (see equation (|l~8|)) for a = | 
versus |m — n|: A = 1.9 (<^>) and A = —3.4 (+) , scaled logarithmically. The lines only 
serve as a guide to the eye. 



5. Conclusion 

In the case of a vanishing magnetic field, a = 0, -P(A) reduces to 
P(A) = A. 



If we put this into the equations ( P~5| ) to (|T§), we get the diagonal elements of the 



Green's Function of the Tight Binding Hamiltonian as calculated in [II|. Thus our 
results are consistent in the limit a = 0. 

From another point of view, we could say, that the minimal coupling to the vector 
potential 

p i— > p + eA, 

that is used to obtain the almost-Mathieu equation for electrons in the one-band-model, 
is, within this model, equivalent to the transformation 

A i-> P(A). 
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In a forthcoming paper |I3| the results will be used to solve Dyson's equation for 
a model with isolated impurities, i.e. to calculate impurity states superimposed on 
Hofstadter's butterfly. 
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